Thermal conductivity of a classical one dimensional spin-phonon system 
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We investigate the thermal conductivity k of the classical Heisenberg spin chain coupled to a 
variety of phonon systems using a Green-Kubo approach as well as coupling to heat reservoirs. 
The decay of energy current correlations is power law or exponential, depending on the type of 
(an)harmonicity of the phonon system and the presence of a substrate. In particular, we find that 
coupling the spin system to displacements with an acoustic type harmonic inter-site potential leads 
to a diverging n. Adding a substrate, harmonic or periodic such as sine-Gordon, drives the system to 
diffusive behavior. With anharmonic potential without a substrate the generic transport is diffusive 
even for the integrable Toda potential, while for the FPU chain is ballistic. 



Recent experiments jjj, |2| convincingly promoted mag- 
netic excitations as a very efficient mechanism for energy 
transport in quasi-one dimensional materials. In partic- 
ular, spin-1/2 Heisenberg chain and ladder (undoped or 
hole doped) as well as spin-1 compounds have been stud- 
ied. The observed highly anisotropic thermal conductiv- 
ity, that is comparable in magnitude to that of metallic 
systems, was attributed to magnetic transport. Theoreti- 
cally, it has been noticed that in the one dimensional (ID) 
spin-1/2 Heisenberg model the energy current commutes 
with the Hamiltonian and thus the thermal conductivity 
is ballistic at all temperatures. This observation falls in 
line with a proposal of unconventional transport in ID 
integrable systems|3j. 

In general, the energy transport in magnetic materi- 
als has a phonon (lattice) and a magnetic contribution. 
The thermal conductivity is limited by phonon-phonon 
(anharmonicity) , spin-spin and spin-phonon scattering. 
In the materials discussed above, described by the ID 
spin-1/2 Heisenberg Hamiltonian (assuming that there is 
no disorder), in the absence of spin-spin scattering, the 
thermal conductivity is only limited by phonon-phonon 
and spin-phonon scattering. Clearly there is experimen- 
tal and theoretical interest in understanding the thermal 
conductivity in this singular case[4|, |3, |6(. However, the 
generic case of ID magnetic thermal transport for higher 
or even classical spin (large S limit) is a benchmark that 
is also of experimental interest, e.g. for S = 1,3/2,5/2 
compounds |2j. 

In this work we foucs on the classical Heisenberg model 
coupled to lattice displacements. As both systems are 
classical we can resort to large scale, state of the art nu- 
merical simulations. This is in contrast to the quantum 
problem which is numerically intractable because of the 
unbounded phonon spectrum and thus limited to pertur- 
bative or semi-phenomenological analysis. 

There is a great deal of work done on the thermal con- 
ductivity of ID classical lattices [2j where a variety of be- 
haviors (ballistic vs. diffusive) was established depending 
on the anharmonicity of the lattice and presence of a sub- 



strate. For the classical Heisenberg model it seems that 
the thermal transport is diffusive ja,0|- 

In the following, we present a series of simulations ex- 
ploring the role of the phonon bath transport properties 
to the total thermal conductivity. In particular we find 
that the coupling to acoustic phonons leads to a diverg- 
ing k. This is in contrast to conclusions for the S =1/2 
quantum chain0, Q based on the memory function ap- 
proach. This difference is surprising because we would 
expect that the quantum spin chain (with ballistic trans- 
port) would be more robust to coupling to phonons than 
the classical one. 

We consider the prototype classical Heisenberg spin 
chain Hamiltonian coupled to lattice degrees of freedom 
carrying the spins, 

H = J2 J n,n+l(S n -S n+1 ) + h4(dQ n /dt) 2 

n 

+ Mvl {V([Q n+1 - Q n }/1) + U(Q n /l)} (1) 

where, S„ = Se n (e„ a unit vector), Q n is the displace- 
ment of the nth atom with mass M from the equilibrium 
position nl, I is the lattice spacing, vq is the velocity of 
longitudinal sound, V (U) the dimcnsionless interatomic 
(substrate) potential (V(0) = V'{0) = 0, V"(0) = 1, 
17(0) = U'(0) = 0, 17" (0) = k > 0). The exchange 
constant J n ,n+i depends on the distance between the 
nth and (n + l)th atoms as J n>n +i = J\\ — x(Qn+i — 
Q n )/!], J > with x the dimensionless spin-phonon cou- 
pling. 

The equations of motion for the magnetic and elastic 
subsystems, 

d o o 9H , r d 2 dH 

^ = -S»x— , M- 2 Q n = -—, (2) 

written in dimcnsionless form are, 

e„ = c f [(l - xPn-i)e n -i X e n - (1- xPn)e n X e n+ i],(3) 

U n = CeXifn-1 - fn) + V' (pn) ~ V (p n -l) - U'(u n ), (4) 
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where the dot denotes differentiation with respect to the 
dimensionless time r = vqI/I, p n = u n +i — u n , u n = 
Q n /l is the dimensionless displacement and /„ = (e n • 
e n+1 ). Here we have introduced, the coefficient c t = 
JSI/vq as the ratio of characteristic elastic (l/vo) and 
magnetic (f / JS) times and c e = JS 2 /Mv 2 as the ratio of 
characteristic magnetic (JS 2 ) and elastic (Mvq) energies. 
In the following, we take JS 2 as the unit of energy. 

We are going to consider three types of inter-cite po- 
tentials, (i) the harmonic potential V(p) — p 2 /2, (ii) the 
FPU potential V(p) = p 2 /2 + p 4 /4, (hi) the Toda po- 
tential V(p) = exp(— p) + p — 1 and two type of on-cite 
substrates, (iv) the harmonic U(u) = ku 2 /2 and (v) the 
anharmonic sine-Gordon, U(u) = 1 — cos(27ru). 

The framework for studying the thermal conductivity 
is the linear response Green-Kubo formalism where the 
thermal conductivity k at temperature T is given by the 
energy current J autocorrelation, C(r) = lim CV(r), 

N—i-oo 

Cn(t) = (J(t)J(0))/NT 2 , 

k = Urn / C{s)ds. (5) 
T ^°° Jo 

N is the number of particles in a chain with periodic 
boundary conditions, J(t) = X^=ii«-( r ) i s the total 
energy current, with the local one j n obtained from 
the continuity equation for the energy; j n = {g n -i + 
9n + (un-i + u n )[c e x{e n -i ■ e„) - V'(p n )]}/2 and g n = 
c t c e (l - xPn-i)(l - xPn)(e n -i x e„ • e„ + i). For x = 
we recover the usual expressions for the isolate spin and 
lattice energy current. The averaging (•} is performed 
over thermalized states of the chain. 

To thcrmalize the system we use Langevin dynam- 
ics with the a white Gaussina noise thermostat coupled 
either to the lattice displacements only or to the spin 
and lattice displacements 9] . We verified that both ap- 
proaches give practically the same results for thermody- 
namic (energy, specific heat) as well as dynamic quanti- 
ties. Following thermalization, the integration of equa- 
tions of motion I|3I4JI gives C(r) up to times r ~ 1000 for 
-/V = 4000 spin/atoms in a chain with periodic boundary 
conditions - large enough to practically eliminate finite 
size effects. As C(t) depends significantly on the con- 
crete realization of the thermalized chain, we perform an 
average over about 10 — 10 5 independent initial realiza- 
tions to improve statistics. In all simulations we consider 
the parameters Ct — 1 and c e = 10 that facilitate the 
convergence of the simulations but are also fairly realis- 
tic as the magnetic energy scale is often larger than the 
phonon scale. 

Let's at first discuss the simplest case of a lattice with 
only harmonic inter-site interaction V(p) — p 2 /2 coupled 
to spins. No substrate is present [U(u) — 0) and thus the 
phonons are acoustic with linear low energy dispersion. 
As there is a zero mode with infinite lifetime due to the 
translational symmetry we take care to use thermalized 
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FIG. 1: The power decay of the autocorrelation function 
C(r) ~ r— a in the spin-phonon chain with harmonic inter- 
cite potential (no substrate, U(u) = 0) for \ — 0.2 and differ- 
ent temperatures: T =0.5, 1, 2, 4 (curves 1, 2, 3, 4). In the 
inset, temperature dependence of power a. 

states with zero center of mass motion 5Z„ =1 m„(t) = 0. 

If x was zero, then the thermal conductivity would 
have two independent contributions. The first, from the 
phonon system that is diverging (purely ballistic) as the 
lattice is harmonic and the energy current is a constant 
of motion. The second, from the spins, that is finite (dif- 
fusive) with the spin energy current correlations decaying 
exponentially in time 9]. Thus C(r) would decay expo- 
nentially fast to a constant (due to the phonon part) and 
the total thermal conductivity would of course be diver- 
gent. 

In the spin-phonon coupled system (^/0), as is shown 
in Fig. ^ the simulations indicate that the decay of the 
total energy current autocorrelations becomes power law 
C(r) ~ r~ a with a < 1. Thus the thermal conductiv- 
ity k of the coupled system - diffusive spin plus ballistic 
phonon transport - still diverges. In the inset, the tem- 
perature dependence of a is shown. It is nonmonotonic, 
with a minimum at T ~ 0.5, and a — » 1 as T — > 0, oo, im- 
plying a logarithmically diverging conductivity at these 
limits. 

In order to analyze mechanisms for breaking the ballis- 
tic behavior, we will now study phonon systems with sub- 
strate or with anharmonic inter-site interaction. First, 
we consider the same harmonic inter-site potential but 
now in the presence of a substrate harmonic potential 
U(u) — ku 2 /2 with k = 0.25. Thus the phonon system 
remains harmonic and the corresponding phonon ther- 
mal transport ballistic. The same discussion as above 
applies for the case x — 0. Now however, as the simu- 
lations indicate, again a power law decay of the energy 
current autocorrelations is found C(r) ~ r~ a but with a 
power a > 1 which implies a finite total thermal conduc- 
tivity. The temperature dependence of the spin-phonon 
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FIG. 2: Thermal conductivity k as a function of temperature 
for spin system coupled to harmonic inter-site and substrate 
potential (curve 1). For reference, thermal conductivity K ap 
of uncoupled Heisenberg spin model (curve 2). In the inset, 
temperature dependence of a in C(r) ~ r~ a - straight lines 
are aid to the eye. 



thermal conductivity is shown in Fig. [21 along with the 
thermal conductivity of the spin system n sp alone. We 
note that, unlike the spin thermal conductivity that de- 
creases basically as 1/T 2 at high T due to the bounded 
spin spectrum, the total one increases linearly with T at 
high temperatures (note that the uncoupled phonon part 
diverges) . Thus it seems that the presence of a substrate 
is necessary for obtaining a finite k in the spin-phonon 
coupled system. The temperature dependence of a is 
shown in the inset indicating marginal, logarithmic con- 
vergence of the energy current autocorrelations as T — > 
and a stronger one - albeit still power law - at high T. 
At low temperatures, although the phonon contribution 
is divergent, the total thermal conductivity is lowered 
compared to the spin part. Thus the phonons basically 
act only as scatterers to the spin excitations. Finally, 
note that as c e = 10, the ratio of magnetic to lattice en- 
ergies, the phonon system is highly excited at T / J ~ 1 
(T/Mv 2 ~ 10). 

The above results were found using the Green-Kubo 
approach. Direct simulation of heat transport in a finite 
chain of length N, with ends in contact to heat reser- 
voirs at different temperatures, gives good agreement be- 
tween the two methods. In particular, for the harmonic 
chain without substrate the heat conductivity coefficient 
diverges as k(N) ~ A - 34 for N -> oo and T = 1. For the 
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FIG. 3: Thermal conductivity k as a function of temperature 
for harmonic inter-site and anharmonic (periodic) substrate 
potentials for coupling \ = 0j 0.1, 0.2 (curve 1, 2, 3). The 
case x = corresponds to the sum of uncoupled spin n 3p and 
phonon k p h contributions. In the inset, temperature depen- 
dence of individual K sp and k p h (curve 4 and 5). 



chain with harmonic substrate (k — 0.25) a finite limit 
exists, k = 8.57, compared to the value k = 8.4 obtained 
from (J5J. 

To study the coupling of the diffusive spin system with 
a diffusive phonon one, we consider the on-site sine- 
Gordon potential U(u) = 1 — cos(27ra) with the same 
harmonic inter-site one V(p) = p 2 /2. This anharmonic 
phonon system is known^cj [n| to show normal (diffu- 
sive) thermal conductivity. Thus for the uncoupled sys- 
tem x = 0, the autocorrelations of each component decay 
exponentially with time and the thermal conductivities 
simply add. For x / 0, now the simulations show that 
the energy current autocorrelation of the spin-phonon 
coupled system (as the individual ones) decays exponen- 
tially C(t) ~ e~ 7T , 7 > 0, and the thermal conductivity 
is finite. In the inset of Fig. we show the tempera- 
ture dependence of the individual ones for x = - K sp 
for the spin, K p h for the phonon subsystem. As we see 
Kph increases monotonically with T in the temperature 
region we are considering, while n sp decreases roughly as 
1/T 2 . Actually we expect that at lower temperatures K p h 
increases again as the phonon system becomes asymp- 
totically harmonic [l2^. The thermal conductivity K of 
the coupled system, Fig. |3 goes through a minimum at 
T w 0.5. It is interesting to note that the sum of in- 
dividual component, k for \ = 0, is significantly larger 
than the actual k for \ 7^ indicating the importance of 
spin-phonon scattering in limiting the thermal conduc- 
tivity. The change 6k/ k is largest, about 50 — 80%, at 
low temperatures and decreases saturating to 20 — 30% 
at high temperatures. 

Finally, we will consider anharmonic phonon systems 
without substrate. A classical example is the FPU 
chainQ that has been extensively investigated and shows 
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FIG. 4: The power decay of the autocorrelation function 
C(t) in the spin-phonon chain with FPU inter-cite potential 
(no substrate, U(u) = 0) for T = 0.2 and x = 0.2 (curve 1). 
For comparison, C(r) for the isolated (x = 0) FPU and spin 
chains (curves 2 and 3) also are given. 



infinite conductivity. The autocorrelation function C(r) 
tends to zero as a power law t~" with a < 1 that de- 
pends on temperature T. At high temperature T = 10, 
a = 0.63 [H, while for T = 0.2, a = 0.82 (see Fig. EJ. 
The spin-phonon system turns out to also show ballistic 
behavior. As shown in Fig. a difference to the FPU 
case is that here exist two time scales, t < tq where 
C(t) decreases exponentially fast, turning to a power 
law for t > t . For T = 0.2, the switch time To « 90 
and increases as the temperature tends to zero, a result 
of the dominating spin transport. Note here, that the 
phonon FPU chain with harmonic substrate has finite 
conductivity^) (power law decay of correlations with 
a > 1) and of course the corresponding spin-phonon sys- 
tem has too. 

Another classic anharmonic phonon system is the Toda 
chain, an integrable system with ballistic transport 0|. 
Our simulations (not shown) conclude that the spin- Toda 
chain has finite thermal conductivity with exponentially 
decaying correlations. The temperature dependence of k 
is qualitatively similar to the one shown in Fig. |31 As 
the Toda chain with substrate potential (a nonintegrable 
system) becomes diffusive, we can speculate that the spin 
system acts as a substrate that breaks the integrability. 
Other spin-phonon examples with finite conductivity are 
the anharmonic Morse and Lenard-Jones potentials. 

One surprizing outcome for the composite spin-phonon 
system with a linear vibrational chain is that while con- 
ductivity is infinite when the phonon chain is purely ac- 
coustic, it becomes finite when there is a phonon gap. 
Numerics shows that in the high temperature regime 
thermalized spins act through the term c e x(/n-i — fn) 
of Eq. as dynamic scatterers for phonons with 

a power spectrum similar but not identical to that of 
white noise. Specifically, we find that this effective 



additive noise term has zero mean and is correlated 
with r c = 0.5. In this high-temperature regime one 
may use the analysis of mode coupling theory relating 
the energy autocorrelation to the phonon group veloc- 
ity v g {k) — dw{k)/dk and the effective local damping 
7(k) as C(r) ~ J^dkv 2 (k)e-^ T ^. Algebraic depen- 
dence of damping on the wavevector fc, viz. j(k) ~ k s 
results in algebraic decay in the purely accoustic chain 
as C(r) ~ r- 1 /^, while the presence of the gap results 
in dominant group velocity contribution at some charac- 
teristic wavevector k ~ kg leading asymptotically to an 
exponential decay as C(r) ~ e~ 7 ( fc °) T . This heuristic pic- 
ture is compatible with the numerical evidence presented 
above. 

In conclusion we find that the thermal conductivity of 
a classical Heisenberg spin chain coupled to linear or non- 
linear phonon chains with a substrate is finite. The same 
holds true when the spin system is coupled to anhar- 
monic chains such as the integrable Toda lattice. How- 
ever, when the phonon chains is either purely linear ac- 
coustic, or weakly nonlinear accoustic such as the FPU 
lattice conductivity becomes anomalous. Further the- 
oretical work is needed for the understanding of these 
numerical observations. 
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